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QUALITATIVE PROPERTIES OF THE BALLISTIC TRAJECTORY.* 

By T. H. Gronwaia. 

1. Introduction. In the following, we shall consider the projectile as a 
particle, and refer it to a right-handed s\-stem of rectangular coordinates 
with the origin at the muzzle of the gun, the positive z-axis being hori- 
zontal and directed toward the target and the positive y-axis vertical and 
directed upward. The retardation R due to the resistance of the air is 
introduced in its most general form as a function of the position and 
velocity of the projectile and the time, so t hat R = R(x, y , z, x', y', z', t), 
where x = dxjdt, etc. Denoting by v = \r' 2 + y' 2 + z' 2 the velocity of 
the projectile, the components of the retardation are Rx'/v, Ry'/'v, Rz' v, 
and g being the acceleration of gravity, the differential equations of 
motion are 

Rx' „ Ry' „ Rz' 
x = - — , y" =- T -g, z -. 

with the initial conditions x = y = z = 0, x! = x ' > 0, y' = y ' > 0, 
z' = for t = 0. The first and third equations give (z'/x')' = 0, so that 
z' = const, x', and since z' = z = for t = 0, it follows that z = for 
any t, that is, the trajectory lies in the xy-pl&ne. 
Consequently the velocity is given by 



V = V + y'\ 

and introducing the notation 

7? 

E = E(x, y, x', y', t) = - , 

the differential equations of the trajectory become 

(1) x" = - Ex', y" = -Ey' - g, 

with the initial conditions for t = 

C2) x = y = 0, x' = x ' > 0, y' = y ' > 0. 

For purposes of comparison, let us consider briefly the trajectory in 
vacuum, where E = 0. Equations (1) and (2) then give 

(3) x' = xo', y' = y ' - gt, x = x 't, y = yo't - igt 2 . 
Besides the origin, two points of the trajectory are of special interest: 

* Read before the American Mathematical Society, April 26, 1919 and February 28, 1920. 

44 



QUALITATIVE PROPERTIES OF THE BALLISTIC TRAJECTORY. 45 

the summit, or point of maximum ordinate, and the point of fall, where 
y = 0. The coordinates, velocity components and time will be denoted 
by the subscript s at the summit and by the subscript a at the point of fall. 
By r we denote the angle of slope defined by tan r = dyjdx and 
— x/2 < t < t/2, so that x' = v cos r, y' = v sin r; the value of r for 
t = is called the angle of departure and denoted by a, and for t = t u , 
we write r„ = — u and call « the angle of fall. The abscissa x„ of the 
point of fall is the range of the trajectory, and the corresponding t„ is the 
time of flight. Eliminating t in (3), we obtain 

(4) y = tana-x - -—#, 

and since the summit and the point of fall are defined analytically by 

y,' = and y m = respectively, we have the following well-known 

formulas 

2x 'yo' v 2 sin 2a 

%a, = = , 

9 Q 



-¥-4 



2x„ tan a 



(5) 



9 

a = a, xj = x ', yj = - y ', 



iy</ 2 

y, = 2 ~Y = 56C = igtj = |x„ tan a, 

t 3 = 0, x,' = xo', y» = 0. 

Returning to the case of a resistance different from zero, we shall now 
derive qualitative properties of the trajectory under various hypotheses 
on the resistance. In § 2, hypothesis A, stating essentially that the 
resistance is positive, will be used to obtain a number of inequalities which 
all reduce to equalities of the type of (5) when the resistance is made equal 
to zero.* In § 3, the stronger hypothesis B, which states that the re- 
sistance depends on v and y alone, decreases or is stationary when y 
increases, and increases faster than v when the latter increases, is used to 
obtain additional inequalities of a similar character. 

* Under the special assumption that the resistance is positive and depends on v alone, so that 
E = E(v) > 0, some of these properties have been obtained before. Thus properties I, III, VI 
VII, X, the first inequality in XI, XII, XIV and XV were given by P. de Saint-Robert, Mem. 
Ac. Sc. Torino, ser. 2, vol. 16 (1S55). Properties XIII and the part of XVI stating that the time 
of flight is less on the rising than on the falling branch, were obtained by Zaboudsky in his Exterior 
Ballistics (St. Petersburg, 1S95). 
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Finally, the still stronger hypothesis C, to the effect that the resistance 
is the product of a function of v by a function of y with some further 
restrictions on the mode of increase of these functions, is introduced in § 4 
for the purpose of investigating the existence of maxima and minima of 
the velocity v.* 

2. Positive resistance — Hypothesis A. All inequalities in this paragraph 
will be derived from 

Hypothesis A : For x is 0, x' > 0, t IS and all finite values of y and 
y', the function E = E(x, y, x', y', t) is positive, and the derivatives dE/dx, 
dE/dy, dE/dx' and dE/dy' exist. Moreover, for any positive values of a, b, c 
and d, there exists an M *= M (a, b, c, d) such that the four derivatives and E 
itself are less in absolute value than M for Si x Si a, < x' Si b, 
— c Si y Si c, — d Si y' Si d and (SO but unrestricted upward. 

Since x' > by equation (6) below, it follows from the general existence 
theorem for the solution of a system of differential equations, that the 
system (1) with the initial conditions (2) has then a unique solution which 
is defined (and finite) for any finite positive value of f.f 

From the first of equations (1), we obtain 

. , - r Edt 

x' = x e Jo , 

and since E > 0, it is seen at once that 

I. The horizontal velocity component x' is positive and decreases steadily 
from the value x ' as t increases from zero. 

Since 



(7) x = f x'dt 

Jo 



* In the special case E = E(v), the result that v has at most one minimum but no maximum 
was obtained by de Saint-Robert (I.e.). 

t When E is not bounded, the trajectory may end in a point which is reached in a finite time, 
beyond which x and y are no longer real. To show this, let us transform (1) to the independent 
variable t; we then obtain the following equations, which are derived in all textbooks on ballistics: 

d(v cos t) tPE dt 



dr g ' dr g cos r ' 

dx v 2 dy v 2 tan t 

d~r ~ g ' dr ~ g 

Now make E = gv~ 3 , which becomes infinite for x' = y' = 0; then the first equation is integrable 

immediately and gives 

{v cos t) 2 = (vo cos a) 2 — 2 sin a + 2 sin t. 

Assume (v cos a) 2 < 2 sin a and determine t x by 

sin ti = sin a — J(» cos a) 2 ; 

then v cos t is imaginary for t < n, and the motion ends at the time h determined by 

U = - I \2sin t — sin ti) — j— . 
g J t\ cos 2 t 
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and x' > by I, it follows that 

II. The horizontal distance x increases steadily from zero as t increases 
from zero. 

Eliminating E between the two equations (1), we find 

(8) dt\x') = ~x" 
and integrating, 

fa\ + + y' v°' C' dt r dx 

(9) tan r - tan a = -, - —, g \ -, = - g \ —, 

•C Xq ,/q X ,7o x 

whence, since x' > 0, 

III. The slope tan r decreases as t increases (or as x increases) so that 
the trajectory is concave downward. 

We may replace (9) by 

„^ dy , rdx 

since x' < x ' by I, we have, 

dy , f x dx 1 x 

-T- < tan a - g I — - 2 = tan a - g — 7i 

ax Jo Xo Xq 

and integrating between the limits and x, 

x 2 
(11) y < tana-x - g— 75, 

£x§ 

whence by comparison to (4), 

IV. The trajectory lies below the trajectory in vacuum corresponding to 
the same initial velocity components x ' and y a ' . 

At an extreme of y, we have y' — 0, and consequently y" = — g < 
by (1), so that every extreme is a maximum. Therefore y, being less, 
by IV, than the maximum ordinate yo' 2 j2g of the trajectory in vacuum, 
has a unique maximum and no minimum, so that 

V. The altitude y increases from t = to its maximum y, at t = t 3 , and 
decreases steadily as t increases beyond t s .* 

From the second of equations (1), it is seen that when y' S 0, then 
y" < so that y' decreases as t incre ases, and since x' also decreases, the 
same is the case with v — \'x' 2 + y'~, or 

VI. The velocity and its vertical component decrease as t increases from 
to t s . 

* In case 1/0' < 0, it follows in the same way that y decreases steadily as ( increases from zero. 
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(12) 








\im y - f - 

I—* oo X 




• oo 





or 

VII. The angle of slope t decreases steadily toward — (ir/2) as t increases 
indefinitely,. 

We shall now prove 

VIII. The trajectory intersects the x-axis for one and only one t = t u 
beyond t = t $ , (so that the range x u exists). 

There can be only one point of intersection, the trajectory being 
concave downward (III), and since y > for t = t s , this point of inter- 
section exists if we can show that y -> — « as t -» « . By I, x' decreases 
toward a limit g as ( -* », which leads us to distinguish three cases. 
First, assume that lim x' > 0; then (12) shows that there exists a h 

t — »-0O 

such that y' < — 1 for £ > ^ and consequently y — yi < — (t — ti), or 
y-»— °= as <->«.* Second, assume that lim x' = 0, but that there 

exists an e > and a t\ such that y' < — e when £ > iij then y — y x 
< — e(t — h) and y -» — « as before.f Finally, assume that x' -» and 
y is bounded downward as < -* oo, but that y' approaches zero infinitely 
often, that is, for any e > however small, the set of ^-values defined by 
— « = y' < contains values greater than an arbitrarily chosen t\. 
Since y is bounded, it follows from (11) that x is also bounded, and since 
< x' < x ', hypothesis A shows that there exists an M such that 
\E\ < M for all values of t for which - 1 g t/' < 0. Now take for e 
the smaller of the values 1 and g/2M; from the second equation (1), we 
obtain 

y"S \Ey'\ -g<M-{ M -g= -\g 

for all t such that — e 51 y' < 0. But y' obviously cannot decrease for 
all t in this set; for some point belonging to this set, we must therefore 

* Example: E = x' — a, where < o < x '. The first equation (1) gives at once 

x' — o Xo' — a , 
— , — = ,— e~<", 

X Xo 

whence x' — ► a as ( -*- x . 

t Example: E = 1. Equations (1) give x' = xo'e~', y' + g = (i/o' + g)e~', so that x' — 
and j/' -*■ — 5 as ( — ► » . 
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have y" S 0, in contradiction to y" < — \g. Hence the last case cannot 
occur, and our theorem is proved.* 

In order to derive further inequalities, we shall use the following 
theorem which is due to Tchebychef : 

When u(x) and v(x) are continuous for a Si x ^b, and both these func- 
tions increase (or both decrease) as x increases, then 



but 



nh s\b nib 

(b — a) I u(x)v(x)dx > I u(x)dx- I v(x)dx, 

Jf*»6 pb /*b 

u{x)v{x)dx < I u{x)dx - I v(x)dx 
a */a *Ja 



when one of the functions increases and the other decreases.^ 
Returning to our trajectory, we form the expression 



d_ 
dx 



. , , 1 d ( y'\ x d (y'\ x 

(y -xtanr) = ^ J t [y ~ *-,) = ~^J t {^) = ftp, 



* That some assumption such as Hypothesis A is necessary to prove VIII is seen by the 
example E = 3(1 + 0*> where E increases indefinitely with t. 
From equations (1), 

x' = Xo'e l -( 1+ » s , 



y ' = e-a+o" [ e y ' - g f^ e^^'du ] . 

+ uy < (i+ ty 



When < u < t, we have (1 + u) 3 < (1 + J) 2 (l + w) and consequently 

e U+0 3 

so that 



y' > j/oV-u-*-" 3 - 



V > V' X' *-*■">'* -g(l- —— ) . 

For a sufficiently large t/o', it follows that y > for large values of t, and consequently the tra- 
jectory does not intersect the z-axis (except at the origin). 

t The following very simple proof was given by F. Franklin in the American Journal of 
Mathematics, vol. 7, p. 377 (1884) : Consider the double integral 

XX ^ — "WIW') — vQc)]dtdx = X *X' WW) - "( J Mz) - u(x)v(t) + u{x)v(x)\dx 

Xb /»6 t*b rb pb nb 

u(t)v(t)dt—J u(t)dt-J v(x)dx — J v(t)dt-J u(x)dx + (b — a) J u(x)v(x)dx, 



pb /»* />b nb fb 

I J J [u(t) — u(x)][v(t) — v(x)]dtdx = (b — a) J u(x)v(x)dx — J u(x)dx • J v(x)dx. 



or replacing t by x in the single integrals to the right 
1 C b C b r / /-M,. ,-, C b 

When u{x) and v(x) vary in the same sense, u(t) — u(x) and v(t) — v(x) have the same sign, so 
that their product is positive, and therefore also the double integral, which proves the first part 
of the theorem. The second part follows from the first upon replacing u(x) by — u(x). 
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the last step being a consequence of (8). Integrating from zero to x, we 
find 

C x xdx 
(13) y = x tan r + g I — - 5- . 

Jo a; 

In the integral, x increases and so does 1/x' 2 by I; we may therefore apply 
Tchebychef's theorem and obtain 



If* C x dx 

y > x tan t + g • - • I xdx • ~r 2 



by (9), or reducing, 



= x tan t + g • - • 7T ■ (tan r — tan a) 

" x 2 g K ' 



- > |(tan a + tan t). 

00 



Now move the origin to the point X\, y\ on the trajectory; x, y and a are 
then replaced by x — Xi, y — y\ and n, and our inequality becomes 

(14) I^^Ktann + tanr) 

or 

IX. T/ie arithmetic mean of the slopes at any two points on the trajectory 
is less than the slope of the chord joining these two points.* 

In the remainder of this paragraph, we shall denote by x\, y\ a point 
on the rising branch of the trajectory (where y' > 0), and by x%, y 2 a 
point on the falling branch (where y' < 0), the corresponding times and 
velocity components being t h #/, y\ and U, x 2 ', y 2 respectively. 

* This is also readily proved by the trapezoid formula with remainder term for the evaluation 
of definite integrals: 

f*f(x)dx = ?^p [f(x0 + /&,)] -\J2 (x - xOfo - x)f"{x)dx 

which is verified at once by integrating by parts in the last integral. Since x — Xi and xt — x 
are positive in the interval of integration, it follows that when f"{x) < for Xi < x < Xz, then 



F 2 f(x)dx > ^—^ [/(*,) +f(x 2 )]. 



1 x x J w ' 2 
Now let/(x) = tan r; then f'(x) = - g/x'* by (8), and 

/"(*) = p j (- g/x' 2 ) = 2 ff z"/*' 4 = - 2gElx' 3 

by the first of (1). Since E > 0, x' > 0, we have/"(x) < and the last inequality becomes 

2/2 — 2/1 = I tan t dx > -^5 — - [tan n + tan t 2 ], 

which is identical with (14). 
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Making x = x 2 , y = 2/2 = Vi in (14), we obtain tan n + tan r 2 < 0, or 

(15) - T 2 > Tl, 

and in particular, for 2/1 = 2/2 = 0, 

(16) co > a, 

whence 

X. At two points of equal altitude the slope of the trajectory is numerically 
greater on the falling than on the rising branch ; in particular, the angle of 
fall a> is greater than the angle of departure a. 

Next, make Xi = 2/1 = 0, x = x« and y = y s in (14), whence y,/x, 
> \ tan a, and furthermore Xi = x s , j/i = y s , x = x„ and y = 0, whence 
— yJ(x M — x„) > \ tan o), so that, combining these two inequalities, 

(17) fx a tan a < y s < \(x m — x,) tan u. 

Since tan r = dy/dx, we have 



Xi 



J„ tann' x - * 2 ~ J 



tan t 2 ' 



X ° Xl ~ L tann' * 2 ~ Xs ~ X - 



^, — . 1 c/j, 2 tan t 2 ' 

and for y x = y 2 , by means of (15), 

(18) Xi > X„ — X 2 , X s — Xi ^ x 2 — x s , 

where the equality sign holds only when x x = x s = x 2 . In particular, 
for Xi = x s = x 2 , (18) gives x« > §x w , and combining this with the outer 
terms in (17), we find 

XI. The abscissa x s of the summit of the trajectory satisfies the inequalities 



( 19 ) 2 X w < X s < + „„ ~ _l I ~ ' X„ 



tan 
tan a + tan w 



Introducing x> > |x„ in (17), we find 

(20) |x tan a„ < y s < |x w tan co, 

or 

XII. The maximum ordinate y s lies between the maximum ordinates of the 
two trajectories in vacuum having the same range x„ and the angles of de- 
parture a and a respectively.* 

* For another proof based on Rolle's theorem, see Charbonnier, Balistique exterieure ration- 
nelle, vol. 1 (Paris, Doin, 1907), p. 193. 
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Upon multiplication by 2y', the second of equations (1) may be written 
(y' 2 )' = - 2Ey' 2 - 2gy' or, since dy/dt = y', 

iy' 2 = - 2E y' - ^- 

Integrating from y to y, we find, since y> = 0, 

y' 2 = f"2Ey'dy+2g(y,-y). 

Since y' is positive on the rising and negative on the falling branch, it 
follows that 

y x ' 2 > 2g(y a - y{), y 2 ' 2 < 2g(y a - y 2 ), 

which may also be written, observing that y 2 is negative, 



(21) y x ' > <2g{y a - Vl ), - y 2 ' < <2g{y s - y a ). 

When ?/i = y 2 , (21) gives 

(22) y x ' > - y 2 ' 
or 

XIII. At two points of equal altitude on the trajectory, the vertical ve- 
locity component is numerically greater on the rising than on the falling 
branch. 

Since also x\ > x 2 by I, it follows from (22) that for y\ = y 2 

(23) vi > v 2 , 
or 

XIV. At two points of equal altitude on the trajectory, the velocity is 
greater on the rising than on the falling branch. 

For the arc s of the trajectory, we have ds/dy = 1/sin r; hence, denot- 
ing by s„ the total length of arc, 



31 = ^7V ' s B - s 2 = I — 
j sin n j a — 



dy 2 



sin t 2 

and by (15), for y t = y 2 , 

(24) Si > s„ — s 2 . 

In particular, it follows from (24) that 

XV. The length of arc of the rising branch of the trajectory is greater 
than that of the falling branch. 

On account of dtjdy — ljy', we have 



ts - t 



C'dyi r»- dy 2 

Jyr i/l Jut t/2 
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and by (21) 

ts _ h < r d Vl t2 -t s > r dy * 

J n V2gr(j/ S - j/i) ' } ' s Jy t ^2g(y s - y 2 ) ' 

or performing the integrations 

(25) t.- t ,<^ES, t ,- ti> ^ER. 

In particular, for 2/1 = 2/2 = 



(26) t. < \~ <t u - U, 

or 

XVI. The time of flight on the rising branch is less, and that on the falling 
branch greater, than the corresponding time of flight on a trajectory in vacuum 
with the same maximum ordinate. 

It should be noted in this connection that the empirical inequality 

y» > f sC 

proposed by several ballisticians, cannot be universally true. On the 
contrary, the inequality sign must be reversed when the angle of departure 
is sufficiently small, as will be shown on another occasion, in connection 
with the use of power series in exterior ballistics. 

To obtain inequalities connecting t a with x u , we write 

J"' 2 /" s dt 

dt = x' —, 

and apply Tcheby chef's theorem: 

j r"% C H dt 1 1 

h — h < 7 r x'dt • —, = , T ixi - xi) • - (tan n - tan r 2 ) 

H — ti Jt l Jt l x ii — i\ g 

whence 

(27) h-t x < \J X2 - Xl (tan n - tan t 2 ).* 

Applying (27) to the rising and falling branches of the trajectory, we find 



± ^ jx s tan a l(x a — x,) tan <o 
*. < \ > t a -t s < yj 



' Another proof of (27) is obtained from 



u-n=p d 4 

"*\ X 

by the appUcation of Schwarz' inequality : 

J"*X2 rx% dx 1 

dx . ) — ; = (x 2 — Xi).- (tan ti — tan t 2 ). 
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and by addition 

j x a tan a \ {x a - x s ) tan o> 
(28) t a <yj— j- + V ^ • 

Replacing x 8 by ^a; w increases the expression to the right; in fact, it follows 
from x s > \Xu and tan a < tan w that 

Vx s tan a + V(x w — Xs) tan « — ( Vj x u tan a + V|x u tan a) 

= ( Vx> — V|x„) Vtan a — ( V§£ w — Vx w — x,) Vtan w 
< Vtan a(*Jxs + Vx u — a;, — 2 VJxJ, 

and the last factor is negative by the algebraic identity 

(Vx^ + Vx u - x s + V2xJ(Vx^ + Va^ - x s - V2xJ 

= — (-Jits — ^X a — Xs) 2 . 

Therefore (28) may be replaced by 



On the other hand, we have since x' decreases, 

.*''dt t. 
tan 



C'dt ts 
« = 9 I p<9—,, 

r <co 

X^j Xs I X CtC *\ Xs ^t^j * fs/j 



and multiplying these two inequalities 

t,(t u — t a ) > — —- — " tan a. 

But < s (< u — t s ) < \tj, and x u — x s > x a tan a/ (tan a + tan «) by (19), 
so that 

2 ^ tan 2 a 

4 <° gr(tan « + tan «) ' 
and consequently we find that 

XVII. The time of flight is bounded in terms of the range by the inequalities 

<*» V§ • a LIT^.. < '» < >i ( ^ + ^^ 

3. Resistance depending on t> and y only. Hypothesis B. All proper- 
ties in this paragraph will be derived from the following hypothesis (ex- 
cept XIX, which requires the stronger assumption stated) : 
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Hypothesis B : The function E = E(v, y) depends on v and y alone. 
For v > and all finite values of y it is positive and has derivatives dE/dv 
and dE/dy, and for any positive values of a and c, there exists an M = M(a, c) 
such that the two derivatives and E itself are less in absolute value than M 
for < i> S a, — c S y = c. Moreover, 

BE . dE . 

-T- > 0, »- - Oi 

dv dy ' 

so that E increases when v increases, but does not increase when y increases. 
We shall first derive some inequalities involving the horizontal velocity 
components at the origin, the summit and the point of fall. From the 
first of equations (1), we obtain 

*' - -' - X" **" - X" ^T** 

and similarly 

r»- E(v 2 ,y 2 ) 

X > ~ X « = J ^tan^^ 

For j/i = y 2 , we have tan n < — tan t 2 by X, and i>i > y 2 by XIV, so 
that E(v\, t/i) > £(v2, 2/2) by hypothesis B, and consequently x ' — x,' 
> x,' — xj or 

XVIII. The horizontal velocity component at the summit is less than the 
arithmetic mean of its values at the origin and at the point of fall. 

(30) x 3 ' < l(x ' + xj). 
From the first of equations (1), it follows that 

.rV J x' J„ Jo 2/1 J, i'i sinn' 

and similarly 

x,' f y ' E{v 2 , 1/2) dy 2 

log ; = I : . 

x u Jo t' 2 — sin t 2 

If, for j/i = 2/2 (whence ^1 > v 2 ), we have j&(i'i, t/i)/i>i S E(v 2 , y 2 )lv 2 , it 
follows that log (x '/x/) > log {x,",xj), or 

XIX. When, in addition to hypothesis B, 2?(w, y) is such that E(v, y)[v 
does not decrease when v increases and y remains constant, the horizontal 
velocity component at the summit is less than the geometric mean of its values 
at the origin and at the point of fall: 

(31) x.' < \x7x7". 
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In the equation 

xo' — x,' = I Ex'dt, 

Jo 

x' decreases, and E also decreases by hypothesis B, since v decreases (VI) 
and y increases; hence Tchebychef's theorem gives 

x,' - Xs '>\- f'x'dt- f t 'm=^io g (^ } ), 

or 

(32) Xs < U log xo' - log x7 • 

Similarly the equation 



gives 



or 



r-h-r--**-r E * 

Xs %0 «/o X 'Jo -" 

1 1 1 / Xo_\ tan a 

Xs Xq t s \ X 8 / Q 



(33) ts<xow l08x ;r^ x -^^ 

Xo — x s y 

whence by (32) 

,„ . s / / tan a 

(34) x 3 < x Xs — — 

and since x s > |x u by (19), we have 

(35) *'>«%£ *. 

Further inequalities of this type may be obtained, but they are too com- 
plicated to be of much interest. 

We now proceed to derive some properties of the trajectory when t 
increases indefinitely. 

XX. . For all values of t, the velocity is bounded by 



(36) v := max 



\ Vo > E(v[ 0)} 



Since vv' = x'x" + y'y", we obtain from (1) 

(37) v' = - vE - g V j. 

By VI and XIV, v < v for < t ^ t m and we need therefore only show 
that the upper bound (36) cannot be reached for t > t w . Assume t x to be 
the smallest value of t beyond t u for which the upper bound (36) is 
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reached; since v a < v , it follows that v must either increase through t x or 
have a maximum there, so that »' S at h. But at t h we have v ^ v a 
and y < 0, hence E{v, y) ^ -E^o, 0) by hypothesis B, and 

— y'/v = — sin t < 1 

by VII, so that (37) gives v' < g — vE(v , 0). Since we have assumed 
that v S g/E(v , 0) at t h it follows that v' < for this value of t, and this 
contradiction to v' ^ proves our theorem. 
Equation (8) may be written 

dr g x' 

jT = 7 cos 2 r = — g — , 

dt x y v 2 ' 

and therefore dx/dr = — v 2 /g or 

1 C a 
X = - I # 2 dr. 

Now «> 2 is bounded by XX, and r > — w/2; consequently x is bounded 
and being an increasing function of t, x tends towards the limit given by 
making r = — w/2 in the above integral. Hence 

XXI. The falling branch of the trajectory has the vertical asymptote 



i r a 

x = - I v 2 dr. 
Since x = | x'dr, and x' is positive and decreasing, it follows from the 

Jo 

boundedness of x that 

a;' -* as i -* oo . 

4. Maxima and minima of the velocity. Hypothesis C. It is the main 
purpose of this paragraph to show that, under fairly weak assumptions, 
there exist no extremes of v beyond a certain value of t. If there are in- 
finitely many extremes, the points at which they occur must therefore have 
all their limiting points at finite distance, and by a stronger assumption, 
this possibility may also be excluded, so that v has then only a finite 
number of extremes. 

The first step in this investigation consists in showing that vE -* g as 
t -<- oo . Suppose first that there exist a ti such that for t > t\, vE either 
increases or decreases steadily as t increases. Since vE > 0, it then follows 
that vE either increases indefinitely or tends toward a limit greater than 
or equal to zero. Now — y'/v = sin ( — r) ->■ 1 by VII, and unless 
vE -> g, it then follows from (37) that there exist an e > and a t 2 such 
that either v' > e for t > U, or v' < — e for t > t 2 , according as lim vE < g 
or lim vE > g (including vE -* o° ) . Consequently we have either 
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v > v 2 + e(t — t 2 ) —■ <*> as t -*■ oo, in contradiction to XX, or v < v 2 
— e(t — t 2 ) —■ — o° in contradiction to v ^ 0. Hence vE ->-g under 
the assumption made. 

There remains to be considered the case when extremes of vE occur for 
indefinitely increasing values of t. It is obvious that in this case, lim 
iaivE as t —■ =o equals the inferior limit of the minima of vE, and 
lim sup vE as t — <x> equals the superior limit of the maxima. Now 



— [-(i+i£)("+'*)+ifM 

by (37), so that at an extreme of vE 

Since d#/d?; > 0, dE/dy S by hypothesis £, it follows that for y' S 0, 
the expression to the left is positive and that to the right negative or zero, 
so that the extremes of vE occur for y' < 0. The first factor to the left 
and the right hand member being then both positive, it follows that 
vE + gy'/v > 0, and since dE/dv > 0, we have 



i(«+»*)<-i£--* 



moreover, < — y' < v, so that finally 

(38) -,£ <wt< -,£^»., 

at an extreme of vE. 

On account of — y'jv -* 1 as t -* oo , it follows from (38) that under 
hypothesis B 

(39) lim inf vE ^ g. 

t— »-00 

Before proceeding further,, we shall prove the following 

Lemma. Let y(x) and z(x) be solutions of the differential equations 

dy ., . dz . . 

m>M #ie initial conditions y = z = y for x = x , and assume that f(x, u) 
and <p(x, u) are continuous for x is #o an d \u — y{x) \ < e (or for x g x 
and |m — «(«) | < e), and that moreover fix, u) < <p(x, u) for x S x and 
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u = y{x) {or for x s x and u = z(x)). Then 

y{x) < z(x) 
for x > x . 

Since, by hypothesis, f(x, u) and <p(x, u) are continuous at x = x , 
u = y , and/(x , y ) — <p(x , y ) < 0, it follows that /(a;, y(x)) — <p(x, z{x)) 

< for x = x S x + 5, when 5 is sufficiently small, so that dy/dx < dzjdx 
and consequently y < z for x < x < x v + 8. Now assume that y(x) 

< z(x) for x < x < x h but y(xi) = z(xi). Then the same argument as 
before shows that y(x) < z{x) for Xi < x < Xi + S, where 8 is sufficiently 
small, and writing v(x) = y(x) — z(x), we have for h positive and suffi- 
ciently small, v(xi — h) < 0, fl(#i) = 0, v(xi + h) < 0, so that 

y(si - h) - t'Qi) . n . v(xi + h) — v(xi) 
-h > U > /i 

and for h -0, we find dv/dx = for x = x\. But we have 

= /(xi, y(xi)) - <p(xi, y(xi)) < 

by hypothesis, and this contradiction shows that there is no Xi > x for 
which y(xi) = z(x-i), so that y(x) < z(x) for x > x . 

As an application, write 7? = — y, so that »? — 00 as £ -» °o , and com- 
pare the equation 

obtained from (37) by a change of variable, with the equation 

(41) ^ = |_ W ,_„), 

the initial conditions being v = V = Vi for v = Vi and r?i belonging to a 
point on the falling branch of the trajectory. Since < — y' = t\' < v, 
we have 

l-. v -E(v,- v ) <l-E(v,- v ), 

and by the lemma, it follows that 

(42) v(r,) < V{ri) for 77 > ijl 

The mechanical interpretation of this is obvious, since (41) is the equation 
of motion of a particle projected vertically downward. 

Now suppose that V has a minimum = V 2 for v = V2] then, by (41), 
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VzECVz, — 772) = g. For 772 < y < t] 2 + 5 and 5 sufficiently small, we 
have V > V 2 since V 2 is a minimum, and by hypothesis B, it follows that 
VE(V, - v) > V 2 E(V 2 , -n)m V 2 E(V 2 , - y 2 ) = g, so that, by (41), 
dV/dy < or V decreases, which contradicts the assumption of a mini- 
mum. Hence V has at most one extreme, which is a maximum, and 
therefore, as £ —■ <x> or 77 -» 00 , V tends either to a finite limit or to infinity, 
so that r] sufficiently large, Y = c where c ^ 0. Now assume that 
E(c, y) -*■ 00 as y -* — °o for any positive value of the constant c; then 
if V ^ c > 0, it follows from (41) that dV/dri — — <x> as y ->■ — « and 
consequently F -* — 00 which is impossible. Consequently F -»0 as 
i -»• co and a fortiori v -* by (42), or 

XXII. When in addition to hypothesis B, we assume that E(c, y) -*■ °o 
as y -* — 00 /or any positive value however small of the constant c, then 

v ->■ as t -»- 00 . 

Returning to the proof of vE -► g, we now introduce 
Hypothesis C. The function E(v, y) has the form 

E(v, y) = <?(«,) • H(y), 

where G(v) > and G'{v) exists and is > for v > 0, and is bounded for 
< v < a, where a is as large as we please; moreover H{y) > 0, H(y) —■ °° 
as y -* — 00, H'{y) exists and is negative and bounded for — b < y < b, 
where b is as large as we please, and h(y) = — d [log H(y)]/dy satisfies 
the condition 

e~ cy h(y) -*■ 00 as y —■ — 00 

for any positive constant c however small. Finally, any one of the three 
following assumptions is made: 

1) h{y) is bounded as y -— — 00 ; no additional condition on G(v). 

2) j^y - as y - - 00, and G(0) = 0. 

3) For some constant m S 1, and some positive constant C\, 

h(y) 



[H{y)\ 



2/C»+l) 



as y -»■ — 



00 . 



and G(v) > C\V m for v sufficiently small. 

This form of E(v, y) is the one used in actual computations, G(v) being 
given empirically in form of a table and assumed to satisfy the conditions 
CI. The factor H(y) is introduced to account for the decrease in air 
resistance due to the decrease of the density of the air with increasing 
altitude; the expression formerly used H(y) = (1 — ky) n , where k and 
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n are positive constants, has been replaced recently by H(y) = e~ hv , 
where h is a positive constant. It is seen at once that both these expres- 
sions for H(y) satisfy CI. 

Under hypothesis CI, it follows from (38) and XXII that for the 
maxima of vE 

y' 

vE < - g — + h(y) -v 2 — g as t — » , 

so that lim sup vE ^ g as t —■ » , and together with (39) this gives 
vE -~g. 

Under hypothesis C2, it is sufficient to prove lim sup VE(V, — 77) 2 g 
as 77 -<- 00 , since for the same value of 77 > 771, v < V by (42) and conse- 
quently vE(v, - v) < VE(V, - 77). 

From (41) we obtain 

Since (7(0) = 0, we have (7(7) = (7(7) - (7(0) = 7(7'(07), < $ < 1, 
and since (?' is bounded (7 being bounded, tending toward zero as 77 — - °° ), 
it follows that 

< 0(V) \ VG ' m < k, 

where k is a constant. With the equation of comparison 

(44) H = k\g - &(- V)] 

and the initial conditions £ = 7(7(7) for 77 = 77!, it follows from the 
lemma that 7(7(7) < £ for 77 > 77!, so that all we need to prove is 

(45) lim sup H(— 77) • $ S g. 

The solution of (44) taking the value £1 for 77 = 771 is 

-FkHdr,[ . r- PkHdr, n pkHdr,. "I 

£ = e J i> £1 + I Ae J ii d77 + ^ I A.-e-'ii #77 , 

L Jy\\ Jm J 

where 772 is any constant greater than 77!, and H = H(— 77). Integrating 
by parts in the last integral, we find 

n. PkHdr, n 1 d PkHdr, 

= lH em L + J-n'(=7r* dv ' 
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and consequently 



(46) 



„, TT - fkHdr, ph PkHdr, . 

H(- 17) • | - g = gHe J m J j^e J ^ dr, 

+ He J m £ x + g J ke J ™ dr] - jjTZTZS e m \> 



For rj > rj2, the first term to the right is less than 



C" h P kHdr, , 

J jje J m dr] 



9 

'is 



1 fkHdv 1 F'kHdr,' 

and by the well-known formula 

<p(v) — <p(ii) _ <p'(vi) 

where 773 is some value between 772 and 77, this expression equals 

g H{- V3 ) 
h h(- 773) ' 
H{- 7,3) 

and since 773 > 772, it follows from hypothesis C2 that this expression may 
be made less than any e > however small by taking 772 sufficiently 
large. Having thus fixed 772, the second term to the right in (46) ap- 
proaches zero as 77 -► «, since 

-pkkdr, 
He J ii — as 77 -* «, 

the logarithmic derivative of this expression being 

*<-'>[K-*]- 

which approaches — 00 when 77 — 00 by hypothesis C2. Consequently, 

(46) gives 

lim sup H( - 77) • £ - g ^ « , 

which is equivalent to (45), and hence vE — • g as t — 00. 
Under hypothesis C3, we have 

\ T - G(V) Hi-vXy- c, F- tf (- 77), 
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and comparing V defined by (41), to W defined by 
(47) <gr_X- ClTr - H (-.,) 

and the initial condition W = V = v x for 7? = rj u we have v < V < W 
for v > Vi- But in (47), we have G(W) = d W m so that (7(0) = 0, and 
hypothesis C3 is stronger in respect to H than C2; hence we have, from 
what has been proved under hypothesis C2, 

lim sup CiW m+1 H(- ij) ^ g, 

and consequently, for 77 sufficiently large 

Cl W" +1 H(- n) < 2g. 

Therefore 

^/<m+« h(y) 



^(2/) < PA(y) < w--h(y) < (ff y 



as j/ — — « by C3, and from (38) we conclude that the upper limit of the 
maxima of vE does not exceed g. Hence vE — g under hypothesis C3, 
so that we have the theorem: 

XXIII. Under hypothesis C, the retardation tends toward the limit g on 
the falling branch of the trajectory : 

(48) vE(v, y) -* g as t — « . 
At an extreme of i>, we have v ' = 0, or 

(49) vE = - g V - 
by (37), and 

= - ~[vE + g V -} ■ v> + vEh(y)-y'+ i v (Ey' + g); 
eliminating vE by (49), we find that for v' = 0, 

(50) v" = 9 - i (gx' 2 -v 2 y"h(y)). 

For a minimum of v, we have v" S 0, and since h(y) > by hypothesis C 
and v- > y'-, it follows from (50) that 



(51) %) 7 < 9 



*">(£)' 



or, since v — ~ V> 
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From equations (1), it is readily seen that 

Since Ey' -* - Ev — - g by (48), and ?/' 2 < v 2 is bounded by XX, it 

follows that there exists a positive constant c such that, for ??i sufficiently 

large, 

d , y' 2 c „ 
^lo gz ->2 for „-,„ 

and consequently 



i r> (l r ) 1 < 



L 

Since e CYI h( — rj) —■ °° as 77 -► 00 by hypothesis C, it follows that 

as t -* 00, 



*<*>(£)' 



and comparing this to (51), it is seen that for sufficienlty large values of t, 
no minima of v can occur, and hence no maxima, since these alternate with 
the minima. Therefore 

XXIV. Under hypothesis C, there exists a U such that for t > h, the 
velocity v decreases steadily toward zero as t increases.* 

If v has an infinity of extremes, the points at which they occur must 
consequently have a limiting point ti at finite distance. By Rolle's 
theorem, the zeros of v" separate those of «/, and consequently t\ is a 
limiting point of zeros of v", and similarly of as many of the higher deriva- 
tives as may exist. By continuity, we therefore have v'(ti) = 0, v"(ti) = 0, 
*"'&) =0, •••. 

Assuming the second derivatives of G(v) and H(y) to exist, v'" exists. 
Differentiating (37) twice, substituting the values of x'" and y'" obtained 
from (1), making v' = v" = 0, and eliminating E by v' = or (49), it is 
shown without difficulty that 

WY" = - 2h(y)g^ + h{y)gy' + h( y y g y' z -?£ + ?£- h'(y)g y -^. 

Equating this expression to zero, and making v" = in (50), we obtain 

* Under the assumption (contrary to C) that h(y) ^ 0, it follows from (50) that v" > 0, 
that is, every extreme is a minimum. The velocity has therefore one minimum (since it decreases 
on the rising branch) and no maximum, or else it decreases steadily from t = onward. 
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after some algebraic reductions 

2m > - "f - 0. 

When h'(y) S 0, this equation is impossible, so that XXIV may be 
replaced by the following: 

XXV. When, in addition to hypothesis C, we have h'(y) S 0, the velocity 
v has only a finite number of extremes. 

In particular, this is true under the assumption made in practice, 
h(y) = h = a positive constant. 

Since v decreases on the rising branch, the first extreme is necessarily 
a minimum, and from v — as t — « it follows that the last extreme 
must be a maximum. In many computed trajectories, a minimum 
followed by a maximum occurs for y > 0, or before the projectile strikes the 
ground (the part of the falling branch where y < being of course of no 
practical interest). In no case, however, has a second minimum been 
discovered, and it is therefore plausible enough that v has only one mini- 
mum and one maximum, even for y < 0. The question of proving or 
disproving this conjecture remains open. 
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